Weak lensing effects in the measurement of the dark energy equation of state with LISA 
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The Laser Interferometer Space Antenna's (LISA's) observation of supermassive binary black holes 
(SMBBH) could provide a new tool for precision cosmography. Inclusion of sub-dominant signal 
harmonics in the inspiral signal allows for high-accuracy sky localization, dramatically improving 
the chances of finding the host galaxy and obtaining its redshift. A SMBBH merger can potentially 
have component masses from a wide range (10 5 — 1O 8 M0) over which parameter accuracies vary 
considerably. We perform an in-depth study in order to understand (i) what fraction of possible 
SMBBH mergers allow for sky localization, depending on the parameters of the source, and (ii) 
how accurately w can be measured when the host galaxy can be identified. We also investigate 
how accuracies on all parameters improve when a knowledge of the sky position can be folded into 
the estimation of errors. We find that w can be measured to within a few percent in most cases, 
if the only error in measuring the luminosity distance is due to LISA's instrumental noise and the 
confusion background from Galactic binaries. However, weak lensing-induced errors will severely 
degrade the accuracy with which w can be obtained, emphasizing that methods to mitigate weak 
lensing effects would be required to take advantage of LISA's full potential. 
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I. INTRODUCTION 
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, The Laser Interferometer Space Antenna (LISA) is a 
space-based gravitational- wave detector that can observe the 
merger of supermassive black holes in a binary with signal- 
to-noise ratios (SNRs) of hundreds to thousands. So far, it 
has not been possible to predict the rate of mergers in the 
Universe, and the event rate for LISA, with any precision. 
The initial seed black holes could have masses anywhere in 
the range ~ 100-10 5 M Q depending on the formation sce- 
nario. Collapse of metal-free massive stars at z <; 20 could 
cad to light seeds, with masses of a few hundred solar masses 
l[. On the other hand, gravitational instability of massive 
aroto-galactic disks could lead to the formation at z <; 10 of 
heavy seeds with masses ~ 1O 5 M 0. These black hole seeds 
might grow by prolonged accretion whereby in-falling matter 
has a constant angular momentum direction and spins-up 
the black hole [H, [j]. Alternatively, the seeds might grow 
by chaotic accretion associated with sporadic in-fall of small 
amounts of matter from a fragmented disc [f|. The merger 
history in LISA's past light-cone depends on how black holes 
formed and evolved. If seed black holes were heavy, then one 
expects several tens of mergers per year and LISA will ob- 
serve all of them. If the seeds were small then the rate in 
LISA's past light cone could be two to three times larger, 
but LISA's sensitivity to smaller black holes will be poorer 
and the number of mergers observed is again a few tens per 



year. Of these, a handful of events might be close enough to 
be useful for cosmography Q. 

LISA may open a new era for precision cosmography since 
black hole binaries are self-calibrating standard sirens 
In astronomy, a standard candle is a source whose absolute 
luminosity can be deduced from certain observed properties 
such as the time- variability of its light-curve, spectral charac- 
teristics, etc. Analogously, the term standard siren is used for 
binary black holes, since the way gravitational waves interact 
with matter is more akin to sound waves, although, of course, 
they are transverse waves traveling at the speed of light. 
Black holes are also self- calibrating since they don't need 
other measures of distance to calibrate their luminosity. This 
is because the luminosity of a binary black hole depends only 
on its chirp mass 1 and its luminosity distance from LISA. For 
a chirping binary, i.e., a binary whose gravitational- wave fre- 
quency changes by an observable amount during the course of 
observation, one can deduce the chirp mass and measure its 
amplitude and thereby infer its luminosity distance. There- 
fore, binary black holes can provide a new calibration for 
the high redshift Universe avoiding all the lower rungs of the 
cosmic distance ladder. 

Supermassive binary black hole (SMBBH) mergers are, 
therefore, potential tools for cosmology. Indeed, a single 
SMBBH observation by LISA might already significant! 
constrain the dark energy equation-of-state parameter w 
fill . From the gravitational waveform, one obtains the lumi- 
nosity distance, Z?l- Due to LISA's orbital motion over the 
observation time, one can also obtain an approximate sky po- 
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sition. If this allows for localization of the host galaxy, then 
rcdshift z can be measured. Since the relationship between 
.Dl and z depends sensitively on w (among other cosmolog- 
ical parameters), the latter can then be constrained. One 
of the hurdles in availing of LISA for cosmology was that 
LISA's angular resolution might not be good enough to lo- 
calize the host galaxy - a step that is crucial for obtaining the 
rcdshift of the source. However, more recent work has mit- 
igated this hurdle by showing that the use of the full signal 
waveform, which contains not only the dominant harmonic 
at twice the orbital frequency 2/ or b, but also other harmonics 
fc/orb, k = 1,3,4,..., in parameter estimation can improve 
the angular resolution to a level that enables the localiza- 
tion of the source for a large fraction of systems in LISA's 
band [IcHl4l |. The sub-dominant harmonics, therefore, are 
essential for precision cosmology with LISA. 

The above conclusions have mainly been drawn by com- 
puting the covariancc matrix of the intrinsic and extrinsic 
parameters associated with an SMBBH. A binary consisting 
of non-spinning black holes on a quasi-circular orbit is char- 
acterized by nine parameters: the chirp mass M. and reduced 
mass fj., a radial vector (-Dl, #n, </?n) giving the location of the 
source, the orientation of the angular momentum (0l, at 
a fiducial time, the epoch of merger tc , and the signal's phase 
ipe at that epoch. The information matrix (inverse of the co- 
variancc matrix) tends to be rather ill-conditioned and it is 
necessary to exercise a lot of care in its computation and 
inversion. Several groups have independently confirmed the 
measurement accuracies and it is now widely believed that 
sub-dominant harmonics truly bring about a dramatic im- 
provement in the estimation of parameters 0, [Iol - [l7j . One 
of the implications is that LISA will be able to measure the 
masses of the component black holes very accurately and 
obtain the mass function of (seed) black holes. 

Localizing the source well enough for its host galaxy to 
be identified is crucial for measuring w, since a knowledge 
of the rcdshift is needed. One may then ask how parameter 
estimation in general benefits from a knowledge of the sky 
position. In the full 9-dimcnsional parameter space, #n and 
0n arc partially correlated with the other 7 parameters. The 
associated degeneracies in parameter estimation get broken 
if (#n,</>n) are exactly known. The covariance matrix then 
gets reduced from a9x9toa7x7 matrix, leading to smaller 
uncertainties [Til [l8{ . 

In Ref. El, it was shown how the inclusion of sub- 
dominant signal harmonics allows for a sufficiently good lo- 
calizability of the source in the sky and measurement of the 
luminosity distance that inference of w becomes possible. 
However, in those papers only a limited number of possi- 
ble parameter values were considered. The event rates men- 
tioned above are integrated over a large range of masses, and 
it is important to know how accurately LISA can measure w 
depending on which binaries it observes during its lifetime; 
indeed, the quality of parameter estimation varies widely de- 
pending on the properties of the sources [l2T - [l4| . Here we 
report on an exhaustive study of the relevant part of param- 
eter space. We look at 15 mass pairs with (observed) com- 
ponent masses roughly in the range 10 5 — 10 s M , at a fixed 



rcdshift. For each of these, a Monte-Carlo simulation is per- 
formed with 5000 instances of sky position and orientation 
of the orbital plane. We then compute in what percentage 
of these cases one would be able to estimate sky position 
well enough that host identification should be possible. Us- 
ing only these instances, we calculate the distribution of the 
uncertainties in w. Our Monte-Carlo results are in such a 
form that they allow for easy rescaling to different redshifts 
z; we study what happens at z = 0.55, z = 0.7, and z = 1. 
LISA's instrumental noise will not be the only restriction in 
the measurement of w; weak lensing will affect the determi- 
nation of luminosity distance. To assess weak lensing effects, 
our distance errors are combined in quadrature with a 4% 
additional error due to weak lensing [l9j . 

This paper is structured as follows. In Section [TT] we give a 
brief description of the signal model used in our study. Our 
choice of systems is explained in Section IIIII and we study 
the impact of source localizability on parameter estimation 
in Section IIVI In Section [V] we describe how we go about 
determining the dark energy parameter w, including our cri- 
teria for localizability of the source. In Section I VII we dis- 
cuss the results of our study, giving the fraction of systems 
for which the host galaxy can be localized, and the level at 
which w can be estimated depending on what kind of system 
is observed. We will show that weak lensing-induced errors 
in £>l will severely limit LISA's ability to measure w. We 
conclude in Section fVIII with the message that future studies 
should focus on correcting the effect of weak lensing. 

Throughout this paper we set G = c = 1 unless stated 
otherwise. 



II. SIGNAL MODEL AND LISA CONFIGURATION 

The coalescence of black hole binary systems is commonly 
divided into three successive epochs: inspiral, merger, and 
ringdown. During the inspiral, the two black holes are well 
separated, and the radial inspiral timescale is much larger 
than the orbital timescale. As a consequence, the grav- 
itational wave signal emitted during this regime is well- 
understood analytically, in terms of the post-Newtonian 
(PN) approximation of general relativity. The latter is a 
pcrturbative approach whereby the amplitudes and phases 
of gravitational waveforms are expressed in terms of a char- 
acteristic velocity v (see [2(| and the extensive references 
therein). The dynamics of the binary system will be mod- 
eled very well in this way for many orbits, but eventually 
there comes a point where the PN approximation fails; af- 
ter that a numerical solution to the full Einstein equations 
is called for. This happens when the innermost stable cir- 
cular orbit (ISCO) is reached, after which the black holes 
plunge towards each other to form a single black hole; this is 
referred to as the merger phase of the coalescence. The re- 
sulting black hole then undergoes "ringdown" as it gradually 
settles down to a quiescent Kerr black hole. 

Although most of the signal-to-noise ratio (SNR) is accu- 
mulated during the final stages of inspiral and merger, dis- 
entangling the parameters that characterize the system and 
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extracting physical information relies critically on longer ob- 
servation times, so that it is important to carefully study 
the inspiral phase. In this work we focus on LISA parame- 
ter estimation from observations of just the inspiral process. 
By doing so we can make use of analytical expressions for 
the detector response and its derivatives with respect to the 
different parameters which are needed to perform parame- 
ter estimation. This allows us to implement fast algorithms 
and explore the parameter space in a comprehensive way, 
unlike in numerical simulations where one can only consider 
a single choice of parameters (masses, spins, . . .) at a time. 
In this way we are able to carry out an extensive study of 
LISA's performance, in particular in measuring the dark en- 
ergy equation of state. Nevertheless, by not including the 
merger and ringdown we are missing a fraction of the to- 
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Here x(t) = [2TrMF(t)} 2/3 is the post-Newtonian expansion 
parameter, with F(t) the instantaneous orbital frequency; 
0(x q ) is referred to as qih PN order. For observed compo- 
nent masses mi and m2, M = m\ + m 2 and v = m 1 rri2/M 2 
are, respectively, the observed total mass and the symmetric 
mass ratio, and £>l is the luminosity distance to the source. 
The explicit expressions for ifv ^ U P to 3PN can be found 
in [25j]. We neglect the contribution of spins, so that the 
gravitational waveform can be parametrized by 9 parame- 
ters: luminosity distance, £>l; two angles (</!>n,#n) defining 
the source position; another two angles (#Lj<^l) specifying 
the orientation of the orbital angular momentum; two mass 
parameters; the phase at coalescence, ipc', and the time of 
coalescence, tc- The sky position and orientation angles are 
defined with respect to a solar system barycentric frame, as 
in 

For LISA observation of supermassive black hole binary 
inspirals, most of the SNR accumulates at frequencies be- 
low 10 mHz, in which case it is appropriate to use the low- 
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where fk = f/k, an overdot denotes the derivative with re- 
spect to time, and ipf,k(t(fk)) is given by 

^/,fc(*(/fc)) = 27r/*(/ fc ) - fc #(*(/*)) -kMt(fk)) - tt/4. 

(2.4) 

The waveform is a superposition of harmonics of the orbital 



tal SNR that LISA would observe (especially for the higher 
mass systems), which would improve parameter estimation. 
While there has been recent work on parameter estimation 
with LISA using numerical waveforms that include merger 
and ringdown pH l22l |. some more understanding may be 
needed before embarking on extensive studies (but see the 
semi-analytic approach of |23| ). As far as inspiral itself is 
concerned, Stavridis et al. [24J studied the effect of spin- 
induced precession of the orbital plane on our ability to mea- 
sure w, though without inclusion of higher signal harmonics 
in the waveform. The work presented here is complemen- 
tary, in that it assumes zero spins but does include higher 
harmonics in the analysis. 

The inspiral PN waveforms in the two polarizations h+ 
and h x take the general form 



(2.1) 

I 

frequency approximation to the LISA response function p6| . 
In this approximation the detector can be regarded as two 
independent Michelson interferometers, and the measured 
strain in each of these separately can be written as 

fc«)(t) = ^[F«\t)h + (t)+FW(t)h x (t)] , (2.2) 

where i = I, II labels the two independent interferometers. 
The response functions F^l' and F^' depend on time through 
the sky position and orientation of the source with respect 
to LISA, which vary over the observation time because of 
LISA's orbital motion. The factor \/3/2 is due to the 60° 
angle between the interferometers' arms. 

It is convenient to express the waveform in the Fourier 
domain using the stationary phase approximation [27} . The 
Fourier transform h^ l \f) of the response of detector i then 
takes the form [MEl: 



(2.3) 



frequency (labeled by the index k), and each harmonic has 
PN contributions to the amplitude (labeled by n; note that 
currently one can only go up to n = 6, as no amplitude cor- 
rections are explicitly known beyond 3PN [IBf). As the PN 
order in amplitude is increased, more and more harmonics 
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appear; at 3PN order there arc eight, which is why the in- 
dex fc only runs up to k = 8. Quantities in Eqs. (|2.3|) and 
(|2.4[) with the argument t(fk) denote their values at the time 
when the instantaneous orbital frequency F(t) sweeps past 



the value f/k. A 



(k,n/2[ 



(£) and 
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\k,n/2) 



(t) are the polariza- 



tion amplitudes and phases of the fcth harmonic appearing 
at ?i/2th PN order, ^(t) is the orbital phase of the binary 
and 4>T>(t) is a time-dependent term representing Doppler 
modulation due to LISA's motion around the Sun. Explicit 



expressions for A 



(0 

(k,n/2) 



and 



/2) 



can be found in 



time dependence of these quantities arises through the beam 
pattern functions due to the varying sky position and ori- 
entation of the source relative to the detector [2(| . The ex- 
pression for 0d (t) is given in (26j . For the PN expansions for 
t(f), #(t), F(t), and F(t) we refer to [13]. 

Each harmonic in hu} (/) is taken to be zero outside a 
certain frequency range. The upper cut-off frequencies are 
dictated by the ISCO, beyond which the PN approximation 
breaks down. For simplicity we assume that this occurs when 
the orbital frequency F(t) reaches -Pisco, the orbital fre- 
quency at ISCO of a test particle in Schwarzschild geome- 
try 2 : Fisco = (6 3 / 2 27rAf Consequently, in the frequency 
domain, the contribution to (J) from the fcth harmonic 
is set to zero for frequencies above fcFisco- Thus, the fcth 
harmonic ends at a frequency 



F 



(k) _ 



ISCO 



2.198 x 10 _3 fc 



10 b Mi 







M 



Hz. 



(2.5) 



In determining the lower cut-off frequencies we assume that 
the source is observed for at most one year, and the fcth 
harmonic is truncated below a frequency kF ln , where Fj n is 
the value of the orbital frequency one year before ISCO is 
reached 



-Fin — F(tisco — Af bs) — 



F 



ISCO 



r + 5M ^"obs^isCO 



3/8 ' 



(2.6) 



For simplicity, Fin was computed using the quadrupole for- 
mula. In the above, tisco an d «isco = l/Vo arc, respec- 
tively, the time and orbital velocity at the last stable or- 
bit, and At b s = 1 yr. However, LISA's sensitivity becomes 
poorer and poorer as one goes to lower frequencies, and cur- 
rent estimates normally assume a "noise wall" at a frequency 
no lower than / s = 10~ 5 Hz. Thus, we take the lower cut-off 
frequency of the fcth harmonic to be the maximum of / s and 

fcF in . 

As has been shown by several gr oups, b oth for Earth-based 
detectors [HEH and for LISA [MIlMllIlt taking into ac- 
count all the harmonics significantly improves the parameter 
estimation, and at the same time it extends the mass reach 



2 Note that the cut-off is placed on the orbital frequency of the binary, 
not on the dominant harmonic in the gravitational wave signal, hence 
the extra factor of 2 in the denominator of the expression for i*isco 
compared to what one often finds in other literature. 



to higher mass systems. However, in our computer code we 
have restricted ourselves to 2PN order in both amplitude and 
phase. We emphasize that there is no technical difficulty in 
going to higher orders, but as shown in fToMTsl [29j . the main 
improvement in parameter estimation occurs in going from 
OPN to 0.5PN order in amplitude, and 2PN order will be 
more than sufficient for our purposes. 

Due to the large SNR values that will be measured by LISA 
in observing SMBBH events, the Fisher information matrix 
formalism can be used to perform the parameter estimation. 
In the limit of high SNR, the probability density distribu- 
tion of the true parameters near the measured value can be 
approximated by a multivariate Gaussian distribution whose 
covariance matrix is given by the inverse of the Fisher infor- 
mation matrix I32| . For each of the interferometers i = I, II, 



the Fisher matrix takes the form 



-.(«) 

a/9 



= 2 



{d a h®\dph®) 



Sn(f) 



df, 



(2.7) 



where d a = d/d\ a , with A Q the parameters to be estimated. 
Specifically, we take these to be 

A = (< c ,0c,cos(6'N),0N,cos(6' L ),(?!> L ,lnL» L , M,fi). (2.8) 

The Fisher matrix for LISA as a whole is then 



r a/ 3 — r^g + r^;. 



(2.9) 



The covariance matrix is E Q/3 = (r 1 ) a ' 9 , which gives the 
covariances between parameters, 



(SX a 5X p ) = S Q/3 , 



and hence also the 1-c uncertainties, 



AA" = [((S\ a ) 2 )] 



1/2 



(2.10) 



(2.11) 



where no summation over repeated indices is assumed. 

As mentioned above, we use the stationary phase approx- 
imation to the gravitational waveform in the frequency do- 
main, which provides a way of getting analytical expressions 
for the observed signals, h^'(f), and for their derivatives 
with respect to the parameters, d a h^(f) [HI [Hi- Only the 
final integral over / in Eq. (|2.7p needs to be done numerically. 

The code we use to generate our results has been vali- 
dated by the LISA Parameter Estimation (LISA PE) Task- 
force [33[ through cross-checking of the output with different 
codes from other groups @. 

In Eq. (|2.7j) . S n (f) is the one-sided noise power spectral 
density, which is a combination of instrumental and galactic 
confusion noise. We take our noise curve to be the one that 
was used by all the members of the LISA PE Taskforce in [f| , 
which also corresponds to the noise curve from the second 
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round of the Mock LISA Data Challenges [3J, |35|. The sky- 
averaged instrumental noise is defined by 



>5'inst(/) — -j~2 



l+ l -( L 

2 V/* 



O.f mHzJ 



f 



4Sg 



(2.12) 



m 

-22 



5 x 10 9 m is the armlength, 
Hz -1 is the (white) position noise level, 
4 Hz -1 is the white acceleration noise 



Hz, L = 

™2 TT„-1 



where / is 
S p = 4 x 10 
S a = 9 x 10 

level and f* = c/ (2irL) is the LISA arm transfer frequency 
(see Ref. [ff for further comments and details). 

The galactic confusion noise is estimated by simulating the 
population synthesis of galactic binaries with periods shorter 
than 2 x 10~ 4 s [36l The confusion noise at the output 
can be fitted as 13811 



<Sconf(/) 



\ f < 10~ 3 Hz 

10" 3 Hz < f < 10" 2 7 Hz 
10- 2 ' 7 Hz < / < 10- 2 4 Hz 
10- 2 ' 4 Hz < / < 10- 2 Hz 
/ > 10~ 2 Hz 

(2.13) 

The total noise curve is the sum of instrumental and confu- 
sion noise, 
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S n (f) = S inst (f) + S coni (f). 



(2.14) 



Finally, as in [6[ , we also apply a lower frequency cut-off at 
10" 5 Hz. 



III. CHOICE OF SYSTEMS STUDIED 

LISA's sensitivity band stretches from 10~ 5 Hz to 0.1 Hz, 
which means (see Eq. (|2.5[) ) that LISA will be able to see the 
coalescence of SMBBH systems of (observed) total mass from 
~ 10 5 Mq to ~ 10 8 Mq with signal-to-noise ratios of several 
hundreds to thousands (sec Fig. [1]) almost anywhere in the 
observable Universe. Signals from higher mass systems will 
not significantly enter the frequency band, and systems with 
a total mass lower than ~ 10 5 Mq will have signal amplitudes 
in LISA's band that quickly become too low to be observable. 

How many SMBBH events is LISA expected to see in a 
year, and of what kind? There are several possible SMBBH 
formation scenarios that are able to reproduce the measured 
optical luminosity function of Active Galactic Nuclei in the 
rcdshift range 1 < z < 6, but they differ in (i) the formation 
mechanism and masses of the "seed" black holes, as well as in 
(ii) the details of how accretion causes black holes to grow in 
time. For instance, Volonteri et al. [l|, [39|, [4(| consider a sce- 
nario where light "seed" black holes (of a few hundred Mq) 
were produced as remnants of metal-free stars at z > 20. Al- 
ternatively, gravitational instability of massive proto-galactic 
disks at z > 10 could have led to the formation of much 
heavier seeds with masses ~ 1O 5 M 0. Regardless of seed 



masses, the seeds may have grown by an accretion process 
in which infalling matter has a constant angular momentum 
direction, spinning up the black holes or by chaotic ac- 

cretion from a fragmented disc, causing much smaller spins 

The implications for LISA of these various scenarios were 
recently assessed by the LISA PE Taskforce [(|. Generally, 
scenarios with heavy seeds lead to several tens of mergers per 
year in LISA's past light cone, all of which will be detectable. 
If the seeds were smaller then the rate in the past light cone 
will be a factor of several larger, but LISA's sensitivity to 
light black hole coalescences is smaller; here too the num- 
ber of observed mergers is a few tens per year. Whatever 
scenario, in the course of its lifetime LISA may see a few 
merger events that are close enough (z < 2) to be useful as 
standard sirens, but their expected masses differ depending 
on the scenario. 

For this reason, in this work we did not focus atten- 
tion on any particular part of the mass range and con- 
sidered a uniform sampling of systems within LISA's mass 
reach. In particular, we have analyzed 15 different pairs of 
masses forming an almost uniform grid in the log mi — logrri2 
plane that covers most 3 of the region of systems observ- 
able by LISA, as shown in Fig. [1] we consider all pos- 
sible combinations of observed component masses mi^ = 
{3.6 x 10 5 , 1.2 x 10 6 , 3.6 x 10 6 , 1.2 x 10 7 , 3.6 x 10 7 }Mq. As 
pointed out by a number of authors, e.g. [1(1 [l2|, 1 3 1 j ] . the 
observed SNRs and parameter estimation depend sensitively 
on the sky location and orientation of the source, so that 
Monte-Carlo (MC) simulations are called for if one wants to 
draw general conclusions about LISA's performance. Thus, 
in this work for each of our 15 pairs of masses, we have car- 
ried out 5000 MCs over {cos#n, </>n, cos #l, </>l} drawn from 
a uniform distribution, computing in every case SNRs and 
parameter errors. 

The luminosity distance, Z?l, and redshift, z are extrinsic 
parameters in the problem. Since SNR, parameter uncer- 
tainties, and masses scale with these in simple ways, one can 
obtain a variety of results from calculations that were done 
with fixed values of Dl and z. 4 

Let us discuss the dependence of the results (SNRs and 
parameter errors) on the luminosity distance and redshift. 
First note that: 

1. From conservation of energy, the amplitude of the grav- 
itational wave signal is inversely proportional to Dl- 

2. The frequency is redshifted because of the expansion 
of the Universe, which causes a blue shift in the ob- 
served masses relative to the physical ones: m b s = 



3 The lower mass systems, despite their low amplitude, remain in the 
LISA band for many cycles, which makes the parameter estimation 
highly expensive in terms of computational time. For this reason we 
restricted ourselves to systems with masses higher than 3.6 X 10 5 Mq. 

4 Note that the output of a simulation is obtained without making any 
assumptions about the relationship between z and Dl- 
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FIG. 1: Left: Median values of the measured signal-to-noise ratio (SNRig P c) from the inspiral phase of SMBBH systems at 1 Gpc, 
as a function of the observed masses of the individual black holes, obtained from 1000 Monte-Carlo simulations over the sky location 
and orientation of the source. The quantity plotted is related to the observed SNR at any luminosity distance Dl as SNRd l = 
SNRi Gpc x^. The superimposed diamonds represent the grid of observed mass cases considered to study parameter estimation and 
its implications in measuring the dark energy equation of state. Right: Representation of how the grid of observed masses in Fig. [T] 
(here as dots) translates into physical masses as we increase the redshift, z, of the source. 



(1 + z)m p hy S , where m p h ys stands for any intrinsic pa- 
rameter with dimensions of mass. 



IV. IMPACT OF SOURCE LOCALIZ ABILITY ON 
PARAMETER ESTIMATION 



Since the SNR and the Fisher information matrix are pro- 
portional to the signal amplitude and its square, respectively, 
we have 5 



SNRoc 



Di 



AD L cx Dl; AA oc D L , 



(3.1) 



where A is any parameter different from Dl- Thus, the quan- 
tities SNR x Dl, ADl/Dl, and AA/D L are independent of 
D L - 

On the other hand, the redshift experienced by any signal 
due to the expansion of the Universe translates into a shift in 
the physical masses. Thus, a single simulation made for some 
pair of observed masses, say {mi.obs , m-2,obs}, will be rep- 
resentative of an infinite set of systems at different redshifts, 
z, with physical masses 



^1, phys 



1 + z ■ 



TT12, phys 



W2,obs 
l + Z ' 



(3.2) 



In Fig. [TJ the panel on the right illustrates how the observed 
masses in the left panel correspond to progressively lower 
physical masses as we consider sources at successively higher 
redshifts. 



5 The error in the luminosity distance behaves differently because we 
are using the derivative with respect to the parameter we are varying. 



Much of the literature on parameter estimation with LISA 
has used the restricted post- Newtonian waveform, which sug- 
gested that the position uncertainty would be too large and 
it would not be possible to find the host galaxy. But as 
pointed out by Sintes and Vecchio [HI, [l6j] and Hellings and 
Moore 1171 . and recently studied more thoroughly by Arun 
et al. [lOjl , Trias and Sintes [H], Ha ], and Porter and Cor- 
nish [14| . the higher harmonics in the orbital frequency that 
will also be present in the signal, carry a significant amount 
of information, and including them in search templates can 
vastly improve parameter estimation. (This is also the case 
for ground-based detectors; see In particular, with the 

inclusion of sub-dominant signal harmonics, in many cases 
the uncertainty in sky position decreases dramatically, so 
that host identification becomes possible, allowing for accu- 
rate measurement of w, as shown in [lfl However, in 
the latter papers only a small number of example systems 
were considered. In the present work we aim to sample the 
parameter space far more thoroughly, enabling a much more 
detailed assessment of what might be possible. 

If the host galaxy of an inspiral event can be found, then 
the sky position will be known with essentially no error. One 
would then match-filter the signal against a template family 
with fixed values for #n and </>n- This will help in remov- 
ing the correlation between (#n,^n) and Dl, resulting in a 
smaller uncertainty in the estimation of the luminosity dis- 
tance than before [ll|. The distance error ADl resulting 
from this smaller Fisher matrix is what determines the error 
on the dark energy equation-of-state parameter w [H| . 

The source localizability criterion that we consider in this 
paper is explained in next section, but first we are interested 
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FIG. 2: Examples of distribution of errors on different parameters before (histograms in the back) and after (front histograms) fixing 
the sky location angles for three different combinations of masses: light and symmetric (left), intermediate and asymmetric (middle), 
and heavy and symmetric (right). AOl, the error on the orientation of the orbital plane, is defined analogously to AJ7n (Eq. I |5.2| l). 
We are assuming a fiducial luminosity distance of Dl = 3 Gpc (zo = 0.55), and we only consider the "localizable" systems for this case. 



Af2 L 

1.533 
3.471 2.179 
4.846 2.947 2.173 
3.581 2.906 2.655 1.722 
14.871 5.937 2.952 2.135 1.279 



1.342 
1.980 1.432 
2.428 1.723 1.507 
1.960 1.753 1.565 1.411 
3.644 2.380 1.750 1.356 1.202 



At c 

1.000 
1.019 1.138 
1.495 1.038 1.005 
2.963 1.146 1.004 1.002 
15.564 1.982 1.150 1.005 1.005 



A</) C 

1.055 
1.399 1.153 
1.756 1.024 1.004 
1.344 1.026 1.002 1.002 
2.830 1.008 1.010 1.004 1.007 



AA4/A4 

1.000 
1.006 1.099 
1.021 1.031 1.004 
1.029 1.011 1.006 1.002 
1.055 1.015 1.010 1.005 1.006 



1.000 
1.006 1.141 
1.021 1.033 1.004 
1.029 1.009 1.002 1.003 
1.055 1.013 1.003 1.005 1.006 



TABLE I: Improvement factors of the errors on different parameters after fixing the sky location angles for the "localizable" systems 
at zo = 0.55. In particular, the quoted numbers represent the ratios between the median values of the error distributions before 
and after fixing the sky location. Note that these are independent of Dl- We provide the results for 6 physical parameters and 15 
choices of observed component mass pairs. The way the values are arranged corresponds to the location of the analyzed systems in the 
log(mi) — log(m2) plane as in Fig. [3] 



in the impact of the localization on the other parameters 
associated with inspiral events. Figure [2] and Table Q] show 
how knowledge of sky position improves the estimation of the 
unit normal to the inspiral plane (where AJIl is the stereal 
angle subtended by the two-dimensional uncertainty ellipse 
on the unit sphere), the luminosity distance Dl, the coales- 
cence time tc, the phase at coalescence <pc, the chirp mass 
A4, and the reduced mass fi. In these figures, we only con- 
sider "localizable" systems, consisting in a certain fraction 
(represented in Fig. [3]) of the total. 
We see the following trends: 

1. Knowing sky position has a bigger effect on light sys- 
tems than on heavy ones. Heavier systems deposit less 



power into LISA's band and have a smaller SNR. In 
that case our localizability requirement can only be 
fulfilled by systems where the correlations between pa- 
rameters were already relatively small to begin with, 
so that parameter estimation would already have been 
good beforehand. Adding the information on sky po- 
sition then will not lead to significantly more improve- 
ment. 

2. Symmetric systems show more improvement in param- 
eter estimation when sky location is known than asym- 
metric ones. Indeed, the odd harmonics are all pro- 
portional to the difference between component masses 
(mi — T112), so that they are absent for symmetric sys- 
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terns. For asymmetric systems, the presence of the 
odd harmonics helps break degeneracies, and adding 
sky position information again does not lead to great 
improvement. 

The improvements in the estimation of chirp mass and re- 
duced mass are modest; depending on the system, the gains 
are between a fraction of a percent and 10% (for M) or 
14% (for fi). Much greater improvements can be seen in 
the measurement of the luminosity distance and the orienta- 
tion of the inspiral plane [i.e., the unit vector determined by 
(#l,0l)], as these are much more strongly correlated with 
sky position. The great accuracy in the determination of Dl 
(typically a fraction of a percent even for quite massive sys- 
tems) will translate into excellent estimation of w, if weak 
lcnsing can be subtracted, as we shall see in the next sections. 



V. COSMOLOGY AND MEASUREMENT OF 
DARK ENERGY 

As argued in |4l[, by treating SMBBH as standard 
sirens, LISA could play a role in studying the physical nature 
of dark energy. From the gravitational wave signal itself one 
can measure the luminosity distance Dl with good accuracy, 
but not the redshift. However, the amplitude and phase mod- 
ulations induced in the observed gravitational waveform due 
to LISA's motion around the Sun allow for a determination 
of the source's position in the sky. If the error ellipse associ- 
ated with the sky position measurement is small enough that 
it contains a sufficiently small number of galaxies or galaxy 
clusters, then it may be feasible to identify the host galaxy, 
possibly with the help of an electromagnetic counterpart to 
the inspiral event. In that case a redshift z can be obtained. 
Now, the relationship between Dl and z depends sensitively 
on cosmological parameters such as Hq, Qm, ^de, and w - 
respectively, the Hubble parameter at the current epoch, the 
matter and dark energy density (normalized by the critical 
density), and the dark energy equation-of-state parameter. 
Hence, separate measurements of the distances and redshifts 
to four or more sources would constrain these parameters. 

For the purposes of this paper we assume a spatially 
flat Fricdman-Lemaitre- Robertson- Walker (FLRW) Universe 
with constant w. In that case, the relationship between the 
luminosity distance Dl and redshift z is given by 

£ L(z)= (i±£) r *l 

#o Jo [n M (i + z'Y + n DE (i + z') 3(1+M,) ] 1/2 

(5.1) 

In principle, a measurement of the cosmological parameters 
could proceed as follows. Imagine that a number of SMBBH 
inspiral events have been found in LISA data, and that their 
host galaxy has been identified. The redshifts z can then 
be determined with negligible error. From the gravitational 
wave signals themselves, the luminosity distances could be 
extracted. A fit of Dl as a function of z using the expression 
(|5.ip would then allow us to deduce the values of Ho, Qmi 
f2rjE 5 and w. In practice, however, there may not be a large 



enough number of sources for which the sky position can 
be determined sufficiently well to allow the identification of 
the host galaxy, in which case it will be impracticable to 
constrain all four cosmological parameters at the same time. 

In this paper we consider an illustrative example. We 
will assume that we only have access to a single inspiral 
event, which will be used to estimate one cosmological pa- 
rameter, w in our case; the other parameters will be con- 
sidered known with negligible errors. Such an observation 
would complement the estimation of (Ho, Qm, ^de, w) ob- 
tained by other means, e.g., via gravitational- wave observa- 
tions of other LISA sources such as the extreme mass ratio 
inspirals (EMRIs) [42[ , or via the observations of stellar mass 
compact binaries with ground-based detectors like the Ad- 
vanced LIGO [H] or the Einstein Telescope [H| . The latter 
may see as many as 500 (stellar mass) inspiral events per 
year with identifiable electromagnetic counterparts, giving 
several thousands over a period of five years. From each 
of the signals a luminosity distance could be extracted, and 
the electromagnetic counterpart would allow us to find the 
host and obtain a value for redshift. Fitting the function 
D^(z) would then severely constrain at least a subset of the 
unknowns (Hq, CIm, ^de, w) [3]. Even without these other 
gravitational wave measurements, by the time LISA is oper- 
ational, all of the parameters (including w) may already have 
been measured with good accuracy through electromagnetic 
means, by continued studies of the Cosmic Microwave Back- 
ground, baryon acoustic oscillations, gravitational lensing, 
and a larger population of Type la supernovae [45[. What 
LISA can add, even if only one parameter is measured, is an 
important consistency check: gravitational wave astronomy 
brings the unique benefit that cosmological parameters can 
be constrained without reference to a cosmic distance ladder. 
It is then natural not to make any a priori assumptions on 
w. 

As explained in the previous section, we simulated a large 
number of instances of SMBBH inspirals, with different 
masses, sky positions, and orientations of the orbital plane 
relative to the observer. We aim to answer two questions: 

1. What fraction of these instances allows for identifica- 
tion of the host galaxy? 

2. For each of the events where the host can be found, 
how accurately can we measure w? 

First we need a criterion to discriminate between cases 
where the host can be identified and cases where it can not. 
In order to do this, we define a fiducial cosmological model 
which we will use as a reference; say, a spatially flat FLRW 
Universe with Ho = 75kms~ 1 Mpc~ 1 , flu = 0.27, Qde = 
0.73, and w = — 1. Using this fiducial model, to the measured 
value of Dl we can associate a fiducial redshift value zq. 
Deciding whether or not the host galaxy can be found will 
involve counting the number of galaxies or galaxy clusters in 
some volume error box around the approximate sky position 
and distance of the inspiral event. The error ellipse in the 
sky, Ar?N, will provide one constraint in determining such 
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FIG. 3: Results for three choices of redshift: z = 0.55 (left column), z = 0.7 (middle column), and z — 1 (right column). In each 
column, the top plots show the fraction of inspiral events for which the host galaxy can be identified (by our criterion N c \ UBtCTS < 3) 
for the various choices of component masses we considered. The bottom plots show the median uncertainties in the determination of 
w for the systems that allow for host determination and hence measurement of redshift. The masses are the physical ones. 



an error box; one has 



2tiV(Acos(6>n)A^n) 2 - (6cos(6» n )<5</>n) 



(5.2) 



However, we are not allowed to also use the error in the 
determination of distance, as measuring w requires an inde- 
pendent determination of the luminosity distance and red- 
shift. Instead, we will count how many galaxies or clusters 
there are in an error box determined by the sky position er- 
ror ellipse AJ7n , and a large redshift interval centered on the 
fiducial redshift zq. For concreteness we take this interval to 



J 



be [0.8zn, 1-2zo]. This is a generous choice: given the above 
values of (Hq, Omj ^de), to reconcile the measured £>l with a 
redshift that differs by 10% from our fiducial zq will typically 
require picking a value of w that lies far outside the existing 
bounds from WMAP and supernovae studies [HI E| • The 
sky position error AJ7n together with the redshift interval 
[0.8zn, 1.2zn] will, through our fiducial cosmological model, 
imply a comoving volume AVc in which to search for host 
galaxies: 



AVc = 



1.220 



dz 



,Aft N Dl(z') 



1 



0.8,0 #o (i + z'f ^/n M (i + Z ')3 + n DE (i + z')^+ w ) ' 

I 



(5.3) 



To arrive at an actual number of galaxies or galaxy clusters of clusters 6 . This density is not known very well; here we fol- 
within the volume AVc we need an estimate for the density 



We note that at redshifts z ~ 1 and beyond, galaxy clusters become 
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low Bahcall et al. (46|, who give p c iusters ~2x 10 -5 /i 3 Mpc~ 3 , 
with h the Hubble parameter at the current epoch in units of 
100 km s~ 1 Mpc~ 1 . The number of clusters within our volume 
error box can then be estimated as -/V c i ust0 rs — /OciustcrsAVc- 
If for a particular inspiral event N c i ustevs turns out to be of 
order 1 then the host cluster can be found and a redshift 
value can be obtained. It could be that the host can be iden- 
tified even when iV c i usters 3> 1, as the binary SMBBH merger 
might be accompanied by a distinctive electromagnetic coun- 
terpart which could be found by future large survey instru- 
ments though electromagnetic counterparts [13, Eg. Even 
so, we will take iV c i us t e rs < 3 to be our localizability crite- 
rion. Issues related to finite cluster size and identification of 
the actual host galaxy will be discussed below. At first in- 
stance one wants to identify the host cluster, and for that it 
will typically not be problematic if there are several clusters 
within the volume box. As an example, consider an inspiral 
at 3 Gpc, which in our fiducial cosmological model corre- 
sponds to z = 0.55. Suppose the volume box contains a few 
clusters with redshifts differing by 10%, e.g., imagine there 
is a potential host cluster at z = 0.6. Then in order to recon- 
cile this slightly larger redshift with the measured distance, 
for the same values of Hq, ^m, and Ode one would arrive at 
w = —0.47, a value that is strongly excluded by WMAP and 
supcrnovae studies. 

When the host galaxy or galaxy cluster of an inspiral event 
can be found, we can assume that the sky position will be 
known with essentially no error, 7 which will allow us to re- 
compute a reduced Fisher matrix removing the correlations 
between (#n ; 0n) and the other parameters. In particular, 
this will translate into an improvement in the estimation of 
.Dl (see Table [TJ) , which is what determines the error on the 
equation-of-state parameter w [Til ]. 

Assuming that (Hq,H,m,^de) are known with sufficient 
accuracy that their uncertainties can be neglected, the error 
on w 0, HH will be given by two contributions: (a) the 
uncertainty in Dl from LISA observations (GW) and (b) 
the error on z from the identification of the host galaxy or 
galaxy cluster (GC), 



Aw = 


dD L 


-l 

D L 




dw 






dD L 


-l 

D L 






dw 





1 dD L 

D L dz 

/ AD L , GC 



-, 1/2 



Az? 



GC 



1/2 



.(5.4) 



The second contribution to Aw could be significant in the 
cases where the host galaxy cluster could be identified, but 



increasingly ill-defined; we will merely use the number of clusters as 
a quantitative way of judging whether the host of an inspiral event 
will be identifiable. 

For systems at z = 1, LISA'S sky resolution will be at best 10 — 4 srad 
. whereas the solid angle subtended by a galaxy cluster at the 
same redshift is at least two orders of magnitude smaller. Hence, 
for the purposes of Fisher matrix calculations, wc can assume that 
the sky position has negligible error when the host cluster can be 
identified. 



not the individual galaxy. Even then, since the typical radius 
of a galaxy cluster is 5 Mpc [49j |. the contribution from this 
term will be 0.17%, 0.12% and 0.08% for z = 0.55, z = 0.7 
and zq = 1, respectively. This is smaller than or comparable 
to the distance error from LISA's noise, so that neglecting it 
would make a difference of at most ~ \f2. In what follows, 
we will assume that the host galaxy can be identified (which 
may be possible with large survey instruments through elec- 
tromagnetic counterparts [13, E3 ) and therefore we neglect 
the contribution from the size of the cluster; but in any case, 
the main conclusions of this article would be the same. 

An important problem in determining distances (also in 
conventional astronomy) is that of weak lensing. As the 
waveform propagates through the matter distribution be- 
tween the source and observer, its amplitude will suffer an 
overall amplification or deamplification, leading to an addi- 
tional random error in the estimation of Dl- For the range 
of distances considered here, this error will be at the level of 
3-5% [jgj . This is substantive: at a distance of 3 Gpc (or 
zq = 0.55), a distance error of 4% will correspond to an un- 
certainty in w of 23%. However, it may be possible to largely 
remove the effect of weak lensing by mapping the mass dis- 
tribution along the line of sight (see, e.g., Refs. [H, HH). In 
the next Section we will see that, if left uncorrected, weak 
lensing will completely dominate over uncertainties due to 
LISA's instrumental noise in the determination of w. 



VI. RESULTS 

We now discuss the results for the localizability of sources, 
and the values of Aw obtained from our Monte-Carlo simu- 
lations. 

The top panels of Fig. [3] show the percentages of "use- 
ful" systems, i.e., the fraction of simulated inspiral events 
for which the sky position error is sufficiently small that the 
host galaxy could be identified, by our criterion iV c i us ters < 3. 
Results are given for three choices of redshift: zq = 0.55, 
zq = 0.7, and z$ = 1, which, in our fiducial model, corre- 
spond to Dl = 3 Gpc, Dl = 4 Gpc and Dl = 6.3 Gpc, 
respectively. Two trends can be seen: 

1 . When the total mass is high, the termination frequency 
2Disco of the dominant harmonic will be low and 
the signal will have less power in LISA's frequency 
band, being "visible" only during the last orbits be- 
fore merger. Both the low observed SNR (see Fig. [T]) 
and the short time over which the signal is observed 
will lead to relatively poor parameter estimation. 

2. Parameter estimation will be worse for symmetric sys- 
tems (i.e., mi = 777-2 ) 5 indeed, all odd harmonics of 
the orbital frequency are proportional to the mass dif- 
ference (mi — 777,2) and hence will vanish in the equal 
mass case. The information they could otherwise have 
carried will then not be present. 

For zo = 0.55, the largest fraction of useful systems is 76.5%, 
which occurs for light and very asymmetric systems with 
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FIG. 4: Accuracies on the estimation of w for Zq = 0.55 = 3 Gpc, according to our fiducial model). The way the plots are arranged 
corresponds to the location of systems in the log(mi) — log(m,2) plane as in Fig. [3] The light (green) distributions are without weak 
lensing; in the dark (red) distributions a 4% error in distance estimation has been folded in, as a heuristic way to account for weak 
lensing. Clearly, if the effects of weak lensing are not removed then they will dominate the uncertainty on w. 



component masses (1.2 x 10 6 ,3.6 x 10 5 )M Q . The smallest 
fraction is 33.5%, for (3.6 x 10 7 ,3.6 x 10 5 )M o ; although 
these systems are even more asymmetric, they are too heavy 
to deposit much power in LISA's band. For zq = 0.7, the 
largest and smallest fractions of useful systems have dropped 
to 62.5% and 11.7%, respectively. For Zq = 1 the analogous 
numbers are 46.0% and 0% respectively. Indeed, for (3.6 x 
10 7 ,3.6 x 1O 7 )M0, there are no systems in our simulated 
population for which A^iusters < 3. 

In the bottom panels of Fig. [31 the median values of 
Au>/|u>| are shown, for the "useful" systems where the host 
can be identified. Binaries with masses (1.2 x 10 6 ,3.6 x 
10 5 )M Q , which gave the largest fraction of useful systems, 
also yield the smallest value for the median error on w: 
for z = 0.55 this is (Aw/|w|) me dian = 0.3%. Still for 



zq = 0.55, the largest median error, (Au>/|u>|) me di a n = 4.1%, 
occurs for very heavy and symmetric systems with masses 
(3.6 x 10 7 ,3.6 x 10 7 )A/ Q . At z = 0.7, the smallest and 
largest median errors on w are 0.4% and 5.5%, respectively. 
For zq = 1, and disregarding the mass pair for which there 
are no useful sources, the smallest and largest median errors 
on w are 0.6% and 8.2%, respectively. 

So far we have discussed errors without taking weak lens- 
ing into account. In Figs. [4j [5l and [U we show the distri- 
butions of errors for each of the mass pairs separately, both 
with and without an additional 4% error folded into the dis- 
tance error through a sum of quadratures, to mimic the effect 
of weak lensing. In all cases, even at zq = 1, LISA's instru- 
mental errors tends to be far smaller than the combined in- 
strumental and weak lensing errors. Indeed, a 4% error in 



12 



1000p 
800 



2000p 
1500 
1000 
500 



-3 -2.5 -2 -1.5 -1 -0.5 
log(A w/|w|) 



2500 




2500 


2000 






2000 


1500 






1500 


1000 






1000 


500 






500 











-2 -1.5 -1 -0.5 

log(Aw/|w|) 



-2 -1.5 -1 -0.5 

log(Aw/|w|) 



2500 
2000 
1500 
1000 
500 





4000 






3000 






2000 






1000 









2500p 
2000 
1500 
1000 
500 



-3 -2.5 -2 



-1.5 -1 -0.5 
lw/|w|) 



-2 -1.5 -1 -0.5 

log(Aw/|w|) 



-2 -1.5 -1 -0.5 

log(Aw/|w|) 



2500p 
2000 
1500 
1000 
500 



2500 
2000 
1500 
1000 
500 



-2 -1.5 -1 -0.5 

log(Aw/|w|) 



-3 -2.5 -2 -1.5 -1 -0.5 

!og(A w/|w|) 



-2 -1.5 -1 -0.5 

log(Aw/|w|) 



-1.5 -1 -0.5 

log(Aw/|w|) 



400 
200 



2000p 
1500 
1000 
500 







1500 






1000 






500 









500p 

400 

300 

200 

100 



_2 -1.5 -1 -0.5 
log(A w/|w|) 




_2 -1.5 -1 -0.5 
log(A w/|w|) 



-3 -2.5 -2 -1.5 -1 -0.5 
log(A w/|w|) 



_2 -1.5 -1 -0.5 
log(A w/|w|) 



-3 -2.5 -2 -1.5 -1 -0.5 

log(A w/|w|) 



FIG. 5: The same as in Fig. H but for z = 0.7 (D L = 4 Gpc). 



£>l translates into an 18.5% to 23% error in w depending on 
redshift, and it is around these high values that the results 
for Aw/|w| are sharply peaked when weak lensing is taken 
into account. This shows that, to make full use of LISA's 
potential, future studies should focus on correcting the weak 
lensing effect. 



VII. CONCLUSIONS 

LISA has the ability to approximately localize supcrmas- 
sive black hole binary coalescence events through modula- 
tion of the observed signal due to LISA's motion around the 
Sun. Electromagnetic follow-up observations could localize 
the source with negligible errors in the source's position in 
the sky and the redshift z of the host galaxy. From the 
gravitational waveform, the luminosity distance can be 
inferred. The relationship between £>l and z depends sen- 



sitively on the past evolution of the Universe, which affects 
the gravitational wave signal as it travels from the source to 
the detector over cosmological distances. Assuming that at 
sufficiently large scales the Universe is approximated well by 
a spatially flat FLRW model and that the Hubble constant, 
the density of matter and the density of dark energy are 
sufficiently well known, the observation of a single SMBBH 
event in LISA can be used to measure the equation-of-state 
parameter of dark energy w. Thus, such events can be used 
as standard sirens, similar to the standard candles of con- 
ventional cosmography, but with no need for calibration of 
distance though a cosmic distance ladder of different kinds 
of sources. 

In [T3, EJ , it was pointed out that inclusion of higher sig- 
nal harmonics has a dramatic effect on source localization, 
making it more likely that we will be able to find the host 
galaxy for relatively close-by {z < 1) SMBBH coalescences; 
these are the ones needed for a good estimation of w. How- 
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FIG. 6: The same as in Fig. H but for z = 1 (D L = 6.3 Gpc). 



ever, in those papers a relatively small number of systems 
were studied. Here we performed large-scale Monte Carlo 
simulations in order to exhaustively probe the relevant part 
of the parameter space. We considered 15 choices of observed 
component mass pairs, each being given 5000 possible sky po- 
sitions and orientations, and the results were scaled to three 
different redshifts. 

For each of the mass and redshift choices, we first com- 
puted how many systems would be localizable by the cri- 
terion that within a generous redshift interval there should 
at most be 3 possible host galaxies or galaxy clusters in the 
sky error ellipse. The fraction of localizable systems varies 
widely with total mass and mass ratio (light and asymmetric 
systems being better), but at our "intermediate" redshift of 
z = 0.7 these were between 11.7% and 62.5%. We note that 
these numbers are likely to be on the conservative side. If 
a coalescence event is accompanied by a sufficiently obvious 



electromagnetic counterpart then our localizability criterion 
may be too strict. Furthermore, a more careful treatment of 
the coalescence process (inclusion of spins as well as merger 
and ringdown) would increase these percentages, as discussed 
below. 

Next we calculated uncertainties on w for those systems 
which passed our localizability requirement so that a redshift 
value would be available. Here too there is considerable de- 
pendence on the mass parameters; at Zq = 0.7 the median 
errors came out to be between 0.4% and 5.5%. 

The waveform model used here included higher signal har- 
monics; indeed, without these it becomes difficult to even 
approximately localize the source in the first place [Tol - [l3| . 
However, spins were ignored, which are also known to im- 
prove parameter estimation [5(| Hl| . Recently Stavridis et 
al. have investigated how well w could be measured with 
the inclusion of spin-induced precession of the orbital plane, 
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though without higher harmonics [24j; they found 1-cr un- 
certainties Aw that are similar to the ones in the present 
work. Although higher harmonics and spins break the same 
degeneracies, no doubt some further improvements can be 
expected by combining the two. We also reiterate that, as in 
Ref. [lU , we only looked at the inspiral signal. Inspiral-only 
position estimates can be improved once again if the merger 



and ringdown signal is taken into account. Babak et al. [21 1 
and Thorpe et al. [22| made a start with this using waveforms 
from numerical relativity simulations, and an extensive study 
on localizability using semi-analytic (non-spinning) inspiral- 
merger-ringdown waveforms was performed by McWilliams 
et al. [23j . It would be interesting to see estimates for 
sky localizability and luminosity distance measurements with 
waveforms that have both higher harmonics and spins in the 
inspiral part, and which incorporate merger and ringdown 
as well. The results concerning localizability of the source 
which we presented here may well be underestimates by fac- 
tors of several. 

As far as the luminosity distance is concerned, if we were 
only limited by LISA's instrumental noise and the confu- 
sion background due to Galactic white dwarf binaries, then 
here too there would be, with no doubt, room for improve- 
ment. However, measurements of Dl get "polluted" by weak 
lcnsing effects, which in turn affects the uncertainties on w. 
Recent work indicates that these effects can be substantially 
reduced by exploiting the brightness of galaxies as a tracer of 



the gravitational fields of matter along the line of sight 52 1 , 
or through mapping shear and flexion of galaxy images 53[ . 
With a deep and wide-field image of galaxies, as might be 
available with, e.g., Extremely Large Telescope [54{ and Eu- 
clid [55l | and with which one could construct a flexion map, 



one may be able to reduce the weak lensing error on £>l 
to about 1.5% in the redshift range 0.5 < z < 1 [56[. At 
z = 0.55 this translates into an 8.5% uncertainty in w, which 
is competitive with electromagnetic measurements. By uti- 
lizing still higher-order effects in the apparent deformation 
of galaxies it may be possible to reduce this number a little 
further. Another idea would be to use high resolution CMB 
maps to estimate weak lensing effects; it is not yet clear, 
however, how one might infer lensing effects at rcdshifts of 
z ~ 1 from such maps. The clear message of the results we 
have presented is that in order to use the full potential of 
LISA as a tool for cosmography, further in-depth studies are 
urgently needed on ways to correct for weak lensing. 



Acknowledgements 

We would like to thank Bernard Schutz for discussions 
at an early stage of this work, and Peter Coles and Mar- 
tin Hendry for discussions on the mitigation of weak lensing 
effects. MT and AMS's work was jointly supported by Euro- 
pean Union FEDER funds, the Spanish Ministry of Science 
and Education (projects FPA2007-60220, HA2007-0042 and 
CSD2009-00064) and by the Govern de les Illes Balears, Con- 
selleria d'Economia, Hisenda i Innovacio. CVDB and BSS 
were partially supported by Science and Technology Facili- 
ties Council, UK, grant PP/F001096/1. CVDB's work was 
also part of the research programme of the Foundation for 
Fundamental Research on Matter (FOM) , which is partially 
supported by the Netherlands Organisation for Scientific Re- 
search (NWO). 



[1] M. Volonteri, F. Haardt and P. Madau, Astrophys. J. 582 
(2003) 559 [arXiv:astro-ph/0207276]. 

[2] M. C. Begelman, M. Volonteri and M. J. Rees, Mon. Not. 
Roy. Astron. Soc. 370 (2006) 289 [arXiv:astro-ph/0602363]. 

[3] J.M. Bardeen, Nature 226 (1970), 64 

[4] K.S. Thorne, Astrophys. J. 191 (1974) 507 

[5] A.R. King and J.E. Pringle, Mon. Not. Roy. As- 
tron. Soc. Lett. 373 (2006) L93 [arXiv:astro-ph/0609598]. 

[6] K. G. Arun et al, Class. Quantum Grav. 26 (2009) 094027 
[arXiv:0811.1011 [gr-qc]]. 

[7] B. F. Schutz, Nature 323, 310 (1986) 

[8] D. E. Holz and S. A. Hughes, Classical and Quant. Gravity, 

20, S65 (2005), [arXiv:astro-ph/0212218]. 
[9] D. E. Holz and S. A. Hughes, Astrophys. J 629, 15 (2005) 

[arXiv: astro-ph/05046 1 6] . 
[10] K. G. Arun et al. Phys. Rev. D 76 (2007) 104016; [Erratum- 

ibid. D 76 (2007) 129903 [arXiv:0707.3920 [astro- ph]]. 
[11] K. G. Arun, C. K. Mishra, C. Van Den Broeck, B. R. Iyer, 

B. S. Sathyaprakash, S. Sinha, Class. Quantum Grav. 26 

(2009) 094021 [arXiv:0810.56727 [gr-qc]] 
[12] M. Trias and A. M. Sintes, Phys. Rev. D 77 (2008) 024030 

[arXiv:0707.4434 [gr-qc]]. 
[13] M. Trias and A. M. Sintes, Class. Quant. Grav. 25 (2008) 

184032 [arXiv:0804.0492 [gr-qc]]. 
[14] E. K. Porter and N. J. Cornish, Phys. Rev. D 78 (2008) 



064005 [arXiv:0804.0332 [gr-qc]] 

[15] A. M. Sintes and A. Vecchio, in Third Amaldi conference 
on Gravitational Waves, edited by S. Meshkov, American 
Institute of Physics Conference Series (American Institute of 
Physics, New York, 2000), p. 403. [arXiv:gr-qc/0005059]. 

[16] A. M. Sintes and A. Vecchio, in Rencontres de 
Moriond: Gravitational Waves and Experimental Grav- 
ity, edited by J. Dumarchez (Frontieres, Paris, 2000). 
[arXiv:gr-qc/0005058]. 

[17] T. A. Moore and R. W. Hellings, Phys. Rev. D 65 (2002) 
062001 [arXiv:gr-qc/9910116]. 

[18] B. F. Schutz, Class. Quantum Grav., 26, 094020 (2009). 

[19] B. Kocsis, Z. Frei, Z. Haiman, and K. Menou, Astro- 
phys. J. 637 (2006) 27 [arXiv:astro-ph/0505394] 

[20] L. Blanchet, Liv. Rev. Rel. 5 (2002) 3 

[21] S. Babak, M. Hannam, S. Husa and B. F. Schutz, 

[arXiv:0806.1591 [gr-qc]]. 
[22] J. I. Thorpe, S. T. McWilliams, B. J. Kelly, R. P. Fahey, 

K. Arnaud and J. G. Baker, Class. Quant. Grav. 26 (2009) 

094026 [arXiv:0811.0833 [astro-ph]]. 
[23] ST. McWilliams, J.I. Thorpe, J.G Baker, and B.J. Kelly 

[arXiv:0911.1078]. 
[24] A. Stavridis, K.G Arun and CM. Will, Phys. Rev. D 80 

(2009) 067501 [arXiv:0907.4686]. 
[25] L. Blanchet, G. Faye, B.R. Iyer, and S. Sinha, Class. Quan- 



15 



turn Grav. 25 (2008) 165003 [arXiv:0802.1249 [gr-qc]]. 
[26] C. Cutler, Phys. Rev. D 57 (1998) 7089 [arXiv:gr- 
qc/9703068]. 

[27] B. S. Sathyaprakash and S. V. Dhurandhar, Phys. Rev. D 44 
(1991) 3819. 

[28] C. Van Den Broeck and A. S. Sengupta, Class. Quant. Grav. 

24 (2007) 155 [arXiv:gr-qc/0607092]. 
[29] C. Van Den Broeck and A. S. Sengupta, Class. Quant. Grav. 

24 (2007) 1089 [arXiv:gr-qc/0610126]. 
[30] L. Blanchet, G. Faye, B. R. Iyer and B. Joguet, Phys. Rev. 

D 65 (2002) 061501; Erratum ibid. D 71 (2005) 129902 

[arXiv:gr-qc/0105099]. 
[31] K. G. Arun et al. Phys. Rev. D 75 (2007) 124002 

[arXiv:0704.1086 [gr-qc]]. 
[32] L. S. Finn, Phys. Rev. D 46, 5236 (1992) [arXivigr- 

qc/9209010]. 

[33] LISA Performance Evaluation wikipage: 

http : //www. tapir . caltech. edu/dokuwiki/lisape :home 
[34] S. Babak et al. [Mock LISA Data Challenge Task Force 

Collaboration], Class. Quant. Grav. 25 (2008) 114037 

[arXiv:0711.2667 [gr-qc]]. 
[35] K. A. Arnaud et al, Class. Quant. Grav. 24 (2007) S551 

[arXiv:gr-qc/0701170]. 
[36] G. Nelemans, L. R. Yungelson and S. F. Portegies Zwart, 

Astron. Astrophys. 375 (2001) 890 [arXiv:astro-ph/0105221]. 
[37] G. Nelemans, L. R. Yungelson and S. F. Portegies Zwart, 

Mon. Not. Roy. Astron. Soc. 349 (2004) 181 [arXiv:astro- 

ph/0312193]. 

[38] N. J. Cornish and T. B. Littenberg, Phys. Rev. D 76 (2007) 
083006 [arXiv:0704.1808 [gr-qc]]. 

[39] A. Sesana, M. Volonteri and F. Haardt, Mon. Not. Roy. As- 
tron. Soc. 377 (2007) 1711 [arXiv:astro-ph/0701556]. 

[40] A. Sesana, M. Volonteri and F. Haardt, Class. Quant. Grav. 
26 (2009) 094033 [arXiv:0810.5554 [astro-ph]]. 



[41] N. Dalai, D. E. Holz, S. A. Hughes, and B. Jain, Phys. Rev. D 

74 (2006) 063006 [arXiv:astro-ph/0601275] 
[42] C. L. MacLeod and C. J. Hogan, Phys. Rev. D 77 (2008) 

043512 [arXiv:0712.0618] 
[43] Samaya Nissanke, Scott A. Hughes, Daniel E. Holz, Neal 

Dalai, Jonathan L. Sievers, [arXiv:0904.1017vl [astro- 

ph.CO]] 

[44] B. F. Schutz, B. S. Sathyaprakash, and C. Van Den Broeck, 

[arXiv:0906.4151 [gr-qc]] 
[45] A. Albrecht et al., [arXiv:astro-ph/0609591] 
[46] N. Bahcall et al., Astrophys. J. Suppl. 148 (2003) 243 
[47] R. N. Lang and S. A. Hughes, Astrophys. J. 667 (2008) 1184 

[arXiv:0710.3795 [astro-ph]] 
[48] B. Kocsis, Z. Haiman, and K. Menou, Astrophys. J. 684 

(2008) 870 

[49] D.W. Hogg, J.G. Cohen and R. Blandford, Astrophys. J. 545 
(2000) 32 

[50] A. Vecchio, Phys. Rev. D 70 (2004) 042001 [arXiv:astro- 
ph/0304051]. 

[51] R.N. Lang and S.A. Hughes, Phys. Rev. D 74 (2006) 122001; 

Erratum ibid. D 77 (2008) 109901 [arXiv:gr-qc/0608062] 
[52] C. Gunnarsson, T. Dahlen, A. Goobar, J. Jonsson, and 

E. Mortsell, Astrophys. J. 640 417 [arXiv:astro-ph/0506764] 
[53] C. Shapiro, D. Bacon, M. Hendry, and B. Hoyle, 

[arXiv:0907.3635vl [astro-ph.CO]] . 
[54] ELT Science Working Group, Report of 

the ELT Science Working Group, 2009; 

http: / / www.eso.org/sci / facilities/eelt /science/doc / 
[55] R. Laureijs et al., ESA Science Document DEM- 

SA-Dc-00001, 2009; http://sci.esa.int/science- 

e /www / object / index. cfm?fobjectid=42822 
[56] M. Hendry, private communication 



